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1.  OBJECTIVES  OF  RESEARCH 

At  UCSD,  we  have  initiated  a  fundamental  research  program  in  order  to  understand 
and  quantify  the  dynamic  response  and  failure  modes  of  rocks  and  concrete,  using  a 
coordinated  effort  which  involves  material  characterization,  high-strain-rate  experiments,  and 
physically-based  analytical-computational  modeling.  The  work  included  major  tasks,  as 
indicated  in  Chart  1. 

Our  work  included  a  set  of  carefully  designed  recovery  experiments  using  UCSD’s 
2.5-  and  6-inch  gas  guns  and  the  1.5-  and  3-inch  Hopkinson  bars;  the  mechanisms  and 
nature  of  shock-induced  damage  in  the  material  are  being  studied.  The  strain  rate,  stress 
amplitude,  and  the  total  input  energy  are  controlled  in  these  tests.  Both  jacketed  and 
unjacketed  samples  are  used.  Through  the  use  of  ultrasonic  measurements,  the  degradation 
in  the  sample  stiffness  will  be  measured  nondestructive^.  These  measurements  will  then  be 
correlated  with  microscopic  observations  of  the  specimens,  using  optical  microscopy  and 
SEM.  Statistical  measures  of  the  microstructural  changes  will  be  established,  using  image- 
processing  techniques.  The  main  objective  at  this  stage  of  the  research  is  to  clearly  identify 
the  microstructure  prior  to  dynamic  loading  and  at  various  stages  of  the  damage 
evolution,  in  order  to  understand  the  physics  of:  (1.)  microcrack  nucleation,  growth, 
interaction,  and  coalescence;  (2.)  plastic  deformation  by  dislocation  motion,  intergranular 
slip,  twinning,  and  void  collapse;  and  (3.)  effects  of  high  strain  rate,  stress  amplitude,  and 
total  energy  input  on  the  nature  and  extent  of  the  induced  damage. 

Based  on  the  knowledge  gained  through  the  above  experimental  observations,  a  set 
of  new  experiments  is  planned  to  study  the  residual  strength  and  dynamic  response  of 
shocked  materials.  Here,  the  sample  geometry  and  the  parameters  of  the  input  shock  will  be 
adjusted,  such  that  shocked  samples  with  different  degrees  of  damage  are  obtained.  These 
experiments  will  be  done,  using  the  6-inch  gas  gun.  The  shocked  samples  will  then  be 
characterized,  and  from  them,  samples  will  be  prepared  for  further  testing.  The  aims  of  this 
second  set  of  experiments  are:  1.  to  quantify  the  residual  strength  and  relate  it  to  the 
preshocked  microstructure  and  the  parameters  which  characterize  the  shock  loading 
conditions;  and  2.  to  identify  the  failure  modes  and  the  failure  micromechanisms  of  the 
preshocked  materials. 

Once  the  dominant  microstructural  features  are  identified,  model  experiments  will  be 
designed  to  directly  examine  the  damage  evolution  process.  Using  transparent  samples 
which  contain  predesigned  defects,  high  strain-rate  microcracking  and  other  related  inelastic 
processes  will  be  recorded  by  high-speed  photography. 

Based  on  our  experimental  observations,  we  are  now  developing  a  physically-based 
micromechanical  to  predict  the  dynamic  response  and  failure  modes  of  both  preshocked  and 
unshocked  materials.  The  predictive  capability  of  these  material-response  models  and  the 
range  of  their  applicability  will  be  tested  by  independent  controlled  dynamic  experiments. 


4/15/00 


2 


UCSD/CEAM 


AFOSR  Final  Technical  Report 


F49620-96- 1-0393 


Chart  1 
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2.  STATUS  OF  RESEARCH  EFFORT .  see  pages  8-32 

3.  PROFESSIONAL  PERSONNEL  ASSOCIATED  WITH  RESEARCH 
EFFORT  AND  DEGREES  AWARDED 

Principal  Investigator:  S.  Nemat-Nasser 

Postdoctoral  Research  Associates:  None 
Visiting  Scholars:  None 

Graduate  Student:  Jacob  Rome 

Degree  Awarded:  Masters  -  Engr  Sciences  (Mechanical  Engineering), 

Spring  1997. 

4.  PUBLICATIONS: 

Work  in  Progress: 

“New  Experimental  Techniques  of  High-Strain  Rate  Testing  of  Concrete” 

Jacob  Rome  and  Sia  Nemat-Nasser 

Abstract: 

We  have  adapted  old  techniques  and  developed  new  techniques  to  study  the  dynamic 
response  of  concrete.  Concrete  was  tested  at  high  strain  rates  using  the  Split 
Pressure  Hopkinson  Bar  (SPHB),  also  known  as  the  Kolsky  Bar.  The  samples 
were  tested  in  both  uniaxial  and  biaxial  compression.  Further,  on  many  uniaxial 
tests,  high  speed  photography  was  employed  to  observe  the  loading  and  unloading 
process  in  the  samples.  Microscopy  and  other  visual  methods  were  employed  to 
observe  the  concrete  after  biaxial  tests. 

5.  INTERACTIONS  (COUPLING  ACTIVITIES) 

5.1  Participation  of  Principal  Investigator  at  Meetings:  Papers 
Presented,  Lectures  and  Seminars 

UCSD’s  IMM  Fourth  Summer  School,  Lecture  presented  “Micromechanics  Approach  to 
Damage”,  July  22  -  August  2, 1996. 

1996  IUTAM  19th  International  Congress  of  Theoretical  and  Applied  Mechanics,  Kyoto 
Japan,  Sectional  Lecture  “Plasticity:  Inelastic  Flow  of  Heterogeneous  Solids  at  Finite 
Strains  and  Rotations,  August  25-31, 1996 

14th  Army  Symposium  on  Solid  Mechanics,  Myrtle  Beach,  SC,  Lecture  presented 
“Fracturing  in  Anisotropic  Brittle  Solids:  Theory  and  Some  Preliminary  Experimental 
Results”,  October  15-17, 1996. 

IMECE’1997,  Atlanta,  GA,  Lectures  Presented  “Role  of  Interface  Phases  on  Failure  of 
Silicon  Nitride”  and  “Crack  Bridging  in  In-Situ  Toughened  Ceramics”,  November  17  -22, 
1996. 
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9th  International  Conference  on  Fracture,  Sydney  Australia,  Invited  Keynote  Lecture 
“Tensile  Fracturing  in  Dynamic  Compression”,  April  1-5, 1997. 

UCSD  IMM  Workshop  on  Multi-Scale  Modeling  of  Polycrystal  Plasticity,  La  Jolla,  CA  , 
Lecture  on  “Dislocation-Based  Models  of  Flow  Stress  of  BCC  and  FCC  Metals  with 
Application  to  Tantalum  and  Copper”,  April  10-11, 1997. 

IMECE’97,  Northwestern  Univ.,  IL,  Organized  IMM  Sponsored  Panel  on  The  Future  of 
Research  Organization  in  the  US:  Universities,  Government  Laboratories,  and  Industry 
Partnerships.  Panel  consisted  of  Drs.  J.  Achenbach  -  Center  for  Quality  Engineering  and 
Failure  Prevention,  NU;  B.  Frost  -  Argonne  National  Laboratory;  and  O.  Richmond  -  Alcoa 
Technology  Center. 

IUTAM  Symposium  on  Micro-  and  Macro-Structure  Aspects  of  Thermoplasticity,  Ruhr- 
Universitat,  Bochum,  Germany,  August  25,  1997,  “Experimentally  Based  Micromechanical 
Modeling  of  Metal  Plasticity  with  Homogenization  from  Micro  -  Macro-scale  Properties”, 
Keynote  Lecture. 

2nd  Euroconference  and  International  Symposium  on  Materials  Instabilities,  University  of 
Thessaloniki,  Greece,  September  1,  1997,  “High- Strain-Rate  Localization  in  Metals: 
Microstructural  Characterization,  Experimental  Observations  and  Computational 
Modeling,”  Invited  Lecture. 

Technical  University  of  Denmark,  Denmark,  September  6,  1997,  “Physically-Based  Metal 
Plasticity:  Rate  and  Temperature  Effects,”  Invited  Lecture. 

ARO  Workshop  on  Mechanics  of  Heterogeneous  Structures,  Research  Triangle  Park,  NC 
March  10-12,  1998,  “Physically-Based  Computational  Modeling  of 

Heterogeneous/Anisotropic  Material  Systems,” 

University  of  Nebraska,  May  29,  1998,  “Dislocation  Barriers  and  Plastic  Flow  of  Metals: 
Experimental  Observation  and  Computation”. 

12th  ASCE  EMD  Conference,  La  Jolla,  CA,  June  17-20,  1998,  “From  Microscale  to 
Continuum  Modeling  of  Frictional  Particulates,”  “Strain  Localization  in  Particulate 
Media,”  “Computational  Modeling  and  Experimental  Plasticity,”  “Bridged  Interface 
Cracks  in  Isotropic  Bi-Materials”. 

13th  National  Congress  of  Applied  Mechanics,  University  of  Florida,  Gainesville,  June  21- 
26,  1998,  “Computational  Modeling  and  Experimental  Plasticity,”  “A  Micromechanical 
Model  for  the  Frictional  Row  of  Particulate  Media”.  Invited  Lecture 

NSF-IMM  Symposium  on  Micromechanical  Modeling  for  Industrial  Materials,  University 
of  Washington,  Seattle,  July  20-22,  1998,  “Multi-Inclusion  Method  for  Finite 
Deformations:  Exact  Results  and  Applications,”  Keynote  Lecture. 

ASME/JSME  Pressure  Vessel  and  Piping  Conference,  San  Diego,  CA,  July  26-30,  1998, 
“Dynamic  Fracture  toughness  of  a  2219-T87  Pressurized  Cylindrical  Vessel,”  “Dynamic 
Fracture  Toughness  of  Miniature  Samples  Using  a  New  Recovery  Hopkinson  Technique”. 

National  Materials  Advisory  Board  -  Guest  Dinner  Speaker,  Sept.  98  “Recent 
Developments  in  Mechanics  of  Advanced  Materials”.  Invited  Lecture 
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IMECE’98,  Anaheim,  CA,  November  1998,  “Localized  Failure  in  Dynamic  Biaxially 
Loaded  Composite”.  Invited  Lecture 

IMECE’98,  Anaheim,  CA,  November  1998,  Presentation  by  Jacob  Rome  (graduate  student 
researcher)  “Concrete  Deformation  under  High  Strain  Rates:  Experimental  Methods  and 
Results,”  Symposium  on  Material  Models  for  Large  Deformations. 

South  Dakota  School  of  Mines  and  Technology,  December  1998,  “Recent  Developments 
in  Mechanics  of  Advanced  Materials”.  Invited  Lecture. 

Tuskegee  University,  Tuskegee,  AL,  February,  1999,  “High- Strain-Rate  High-Temperature 
Response  of  Materials”.  Invited  Lecture. 

Arizona  State  University,  February,  1999,  “Improvable  Bounds  on  Overall  Properties  of 
Heterogeneous  Finite  Solids.” Invited  Lecture 

1999  TMS  Annual  Meeting,  San  Diego,  CA,  March,  1999,  “Dynamic  Behavior  of  SiC 
under  Uniaxial  Compression”.  Invited  Lecture 

15th  US  Army  Symposium  on  Solid  Mechanics,  Myrle  Beach,  SC,  “On  Multi-Scale 
Homogenization  and  Average  Field  Theories”. 

Mechanics  of  Heterogeneous  Materials  Symposium,  Grenoble,  France,  June  1999, 
“Anistropy  in  Deformation  of  Granular  Materials”.  Invited  Lecture 

ASME  Joint  AMD/MD  Conference,  Virginia  Polytechnique  Institute  and  State  University, 
June  1999,  “Tri-Axial  Dynamic  Compression  Experiments  Using  a  Novel  Hopkinson 
Technique,”  “Strain  Localization  in  Frictional  Granules,”  and  “High-Temperature,  High- 
Strain-Rate  Receovery  Techniques  Based  on  Kolsky’s  Split  Hopkinson  Bar  Invention”. 


5.2  Consultative  and  Advisory  Functions  to  Other  Agencies, 
Laboratories  and  Universities 

4th  IMM  Young  Investigators  Meeting  on  Critical  Issues  in  Materials  and  Mechanics, 
Boeing  CO,  Seattle,  WA,  Panelist  for  “Materials  Science  in  Mechanics  of  Metals:  Length 
Scale  and  Heterogeneities”,  July  22-August  2, 1996. 

AFOSR  Shock  Physics  Meeting  -  EAFB,  October  18, 1996. 

UCSD  IMM  Workshop  on  Identification  of  Basic  Research  Issues  Arising  from  Industrial 
uses  of  Polymeric  Structural  Materials,  Santa  Barbara,  CA  -  Panelist  Participant,  January 
12-13, 1997. 

AFOSR  Grantee/Contractor  Conference,  Lansdowne,  VA ,  February  2-5, 1997. 

13th  Ceramic  Modeling  Working  Group  Meeting,  APG,  MD,  Organized  through  Institute 
for  Advanced  Technology  (Austin  TX),  Lecture  on  “Confined  Ceramic/W-Penetrator 
Interaction:  Microstructural  Observation  and  Modeling”,  April  24-25, 1997. 

Joint  DOE/DOD/NASA  Pressure  Safety  Workshop,  Lawrence  Livermore  National  Lab, 
December  3, 1997. 
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AFOSR  Tyndall  &  Eglin  AFB  Grantee  Meeting  -  TAFB,  January  17-19, 1998. 

Workshop  on  an  NSF  Initiative  on  Long  Term  Durability  of  Materials  and  Structures: 
Modeling  and  Accelerated  Techniques,  University  of  California,  San  Diego,  February  3, 
1998. 

ARO  Workshop  on  Mechanics  of  Heterogeneous  Structures,  Research  Triangle  Park,  NC 
March  10-12, 1998. 

IMM  Workshop  on  Materials  for  the  Infrastructure,  La  Jolla,  CA,  April  1-3, 1998. 

ATP  Blue  Ribbon  Panel  on  Fracture  Mechanics,  “Fracture  Mechanics  Models  for  Airborne 
Laser’s  Lethality  Program,  Kirkland,  AFB,  April  1998. 

Tuskegee  University  Short  Course  and  Laboratory  Instruction,  June  10-11,  1999.  Short 
Course  entitled  “Dynamic  Experimentation”. 

ATP  Blue  Ribbon  Panel  on  Fracture  Mechanics,  “Fracture  Mechanics  Models  for  Airborne 
Laser’s  Lethality  Program,  Kirkland,  AFB,  September,  1999. 

6.  PROFESSIONAL  HONORS,  AWARDS  AND  FELLOWSHIPS  RECEIVED 
DURING  GRANT  PERIOD 

Sia  Nemat-Nasser  - 

Recipient  of  John  Dove  Isaacs  Chair  in  Natural  Philosophy  1995-2000 
President  of  the  American  Academy  of  Mechanics,  1996-97 
V/Chair  of  the  ASME  Materials  Division,  1996-97 
Technical  University  of  Crete  -  Gold  Medal,  1997 
Chair  of  the  ASME  Materials  Division,  1997-98 
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Status  of  Research  Effort 


1.0  Summary 

In  the  past  years,  we  have  designed  new  experimental  techniques  to  conduct  a  wide  range  of 
experiments,  developed  a  new  micromechanics  based  model,  characterized  the  material  under 
study,  and  implemented  the  model  using  Mathematica  with  the  data  we  have  collected  from  the 
experiments  and  characterization. 

Hopkinson  bars  and  a  20  KIPS  Instron  load  frame  have  been  used  to  conduct  experiments,  on  the 
concrete  and  its  constituent  materials,  at  high  and  low  strain  rates  under  unixial  and  triaxial  com¬ 
pression.  Ultrasonics  have  been  used  for  non-destructive  evaluation  (NDE)  of  the  materials.  Mate¬ 
rial  characterization  has  been  done  using  both  scanning  electron  microscopy  (SEM),  and  standard 
optical  tools.  A  physically  based  model  has  been  developed,  which  is  based  on  principle  of  micro¬ 
mechanics  and  incorporates  damage  parameters.  A  representative  volume  element  (RVE)  is  the 
basis  of  the  model.  In  this  approach,  a  RVE  is  selected  so  it  contains  the  important  material  fea¬ 
tures  that  exist  in  the  material.  The  response  is  then  found  for  the  RVE  based  on  the  properties  of 
each  constituent,  and  the  entire  structure  is  assumed  to  have  the  same  response.  In  the  RVE  for 
one  proposed  model,  the  mortar,  the  aggregates,  the  interfacial  transition  zone  (ITZ)  and  cracks 
are  included.  The  dynamic  effect  and  the  history  are  both  included  via  the  crack  growth.  This 
model  is  implemented  using  the  commercial  software  Mathematica. 

2.0  Materials 

The  G-mix  concrete  (and  its  constituent  materials)  for  these  experiments  were  generously  sup¬ 
plied  by  Tyndall  Air  Force  Base  in  Florida.  Concrete  samples  with  diameters  of  two-  and  three- 
inches  were  used.  In  addition,  Tyndall  AFB  also  supplied  mortar  samples  without  the  largest 
aggregate,  and  samples  of  the  limestone  used  as  the  large  aggregate  in  the  G-mix  concrete.  Figure 
2.1  and  Figure  2.2,  show  examples  of  the  concrete,  mortar  and  limestone  samples. 


Figure  2.1:  Cut-away  view  of  a  concrete  sample  with  a  3”  diameter 
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3.0  Modeling 

The  model  we  have  developed  is  based  upon  the  principles  of  Micromechanics.  A  Micromechan¬ 
ics  based  model  examines  the  physical  structure  of  a  material  to  construct  a  useful  representative 
volume  element  (RVE).  The  proposed  model  incorporates  the  main  features  of  the  concrete.  As 
shown  in  Figure  3.1,  the  limestone  aggregate  (also  referred  to  as  rock)  is  embedded  inside  the 
mortar.  The  modulus  of  each  material  is  noted.  The  mortar  itself  is  inside  of  an  infinitely  extended 
material  that  has  the  modulus,  C,  of  the  overall  structure  (the  concrete).  As  the  material  is  subject 
to  loading,  the  cracks  will  grow  according  to  a  dynamic  crack  growth  model.  By  tracking  the 
length  of  the  cracks  within  the  model  and  relating  that  to  the  strain,  the  response  of  the  structure  is 
found,  and  the  total  damage  to  the  material  can  be  assessed.  The  different  lengths  of  the  cracks 
will  reveal  any  induced  anisotropy. 
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e°  applied 


Figure  3.1:  Sample  RVE  of  concrete.  In  this  case,  two  ellipsoidal  limestone  aggregates  and  three  penny-shaped 
cracks  are  contained  within  the  mortar.  The  mortar  is  contained  within  the  as-yet  unknown  average  properties 
of  the  overall  material  (the  concrete). 

Three  penny-shaped  cracks,  each  parallel  to  a  different  plane,  are  shown  in  the  above  example. 
Elliptical  or  winged  cracks  could  also  be  considered.  As  the  material  is  loaded,  the  cracks 
lengthen  according  to  a  dynamic  crack-growth  law. 

The  modulus  of  the  mortar  and  limestone  have  been  measured  directly  using  ultrasonics  and  dur¬ 
ing  low  strain  rate  loading,  and  the  volume  fraction  of  limestone  has  been  found  using  optical 
techniques.  In  some  cases,  the  effect  of  the  interfacial  transition  zone  (TTZ)  may  be  included  via 
adjustment  of  the  limestone  modulus.  The  number  and  shape  of  limestone  inclusions  in  the  model 
can  be  chosen  based  upon  the  size  and  density  of  the  limestone  aggregates  in  the  concrete.  In  the 
simplest  case,  one  spherical  inclusion  is  used.  In  any  case,  the  inclusions  must  be  elliposidal  to 
satisfy  the  assumptions  in  the  model.  Likewise,  the  mortar  inclusion  (embedded  within  the  overall 
material)  must  be  ellipsoidal,  and  is  chosen  as  spherical  for  now.  The  volume  of  the  inclusions  is 
chosen  such  that  the  volume  fraction  is  consistent  with  that  found  during  the  material  character¬ 
ization.  The  inital  size  and  density  of  the  cracks  are  based  upon  optical  and  experimental  evalua¬ 
tion. 

3.1  Eigenstrains 

To  explain  the  model,  a  few  key  theories  are  reviewed  [Nemat-Nasser,  1999,  Eshelby].  The  con¬ 
cept  of  eigenstrains  was  enabled  by  work  done  by  Eshelby.  An  inclusion  is  embedded  in  a  homog¬ 
enous  domain  (the  matrix),  which  is  subject  to  strains  e°.  If  no  inclusion  were  present,  then  the 
stresses  would  be: 

c°=C:e°,  (3-D 
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where  C  is  the  elastic  modulus.  However,  the  presence  of  the  inclusion  causes  a  disturbance 
strain,  ed(x),  throughout  the  material.  The  new  strains  could  written  as: 

e(x)=e°+ed(x).  (3.2) 

The  stresses  in  the  material  become: 

o(x)=C:e(x)=C:(e°+ed(x))  in  the  matrix  material  (3.3) 

a(x)=CQ:£(x)=Cn:(e0+£d(x))=a°+od(x)  in  the  inclusion  (3.4) 

In  order  to  simplify  the  problem,  an  equivalent  homogenous  solid  can  be  considered.  Here,  the 
solid  is  considered  to  have  only  one  modulus,  C,  which  is  the  modulus  of  matrix  material.  In  order 

for  this  material  to  have  the  proper  stresses  in  the  inclusion  area,  eigenstrains,  £*,  are  introduced  in 
the  inclusion  region: 

a(x)=C:(£(x)-£*(x))=C:(£°+£d(x)-£*(x))=a°+ad(x)  in  the  inclusion  region  (3.5) 

To  provide  a  consistent  solution,  it  is  clear  that  o(x)  from  equations  (3.4)  and  (3.5)  must  be  equiv¬ 
alent,  therefore: 


a°-K?d(x)=C:(£°+£d(x)-£*(x)) 
The  disturbance  strains  can  then  be  written  as: 

£d(x)=S(x;£*(x)) 


(3.6) 


(3.7) 


where  S  is  Eshelby’s  tensor. 

Eshelby  showed  that  when  the  matrix  material  extends  infinitely  and  is  linearly  elastic,  and  the 

inclusion  is  an  ellipsoid,  the  eigenstrains,  £*,  are  constant  in  the  inclusion,  and  S  (Eshelby’s  Ten¬ 
sor)  is  a  function  only  of  the  geometry  of  the  matrix  material  and  the  elastic  parameters  of  the 
matrix.  When  the  matrix  is  isotropic,  it  depends  only  on  the  geometry  and  the  Poisson  ratio.  This 
result,  constant  eigenstrains  and  the  Eshelby  Tensor,  will  be  exploited  in  the  development  of  this 
new  model. 

3.2  Double  inclusion  method 

The  double  inclusion  method  is  a  scheme  that  is  designed  to  find  the  modulus  for  set  of  materials 
and  geometries  [Tanaka;  Nemat-Nasser,  1999].  It  differs  from  methods,  such  as  the  dilute  distri¬ 
bution  or  self-consistent  methods,  by  better  describing  the  interaction  between  an  inclusion  and 
the  surrounding  materiaLThe  dilute  distribution  embeds  in  inclusion  within  the  matrix  material, 
and  thus  ignores  interaction  effects.  The  self-consistent  method  embeds  the  inclusion  in  a  material 
with  the  as-yet  unknown  modulus  of  the  overall  structure.  This  accounts  for  the  interaction 
between  inclusions  in  a  limited  sense. 
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The  double  inclusion  method  differs.  In  the  simplest  case,  one  ellipsoidal  inclusion  with  modulus 

is  embedded  within  an  ellipsoidal  portion  of  the  matrix  material,  CM.  This  outer  ellipsoid  is 
contained  within  an  unbounded  solid  with  elasticity  C.  See  Figure  3.2.  In  the  general  case,  the  dis¬ 
turbance  strains  caused  by  the  inclusion  are  not  constant.  However,  the  average  eigenstrains  can 
be  found  using  Eshelby’s  results  [Tanaka].  Using  this  result,  and  assuming  piecewise  constant 

eigenstrains  (e.g.,  e*1  in  the  matrix,  £*2  in  the  inclusion),  the  overall  modulus  can  be  estimated. 


*  C 

i _ ■ 

Figure  3.2:  Basic  schematic  of  double  inclusion  method 

3.3  Description  of  the  model 

An  initial  estimate  of  the  modulus  for  the  RVE  without  cracks  is  made  first.  The  average  strain  is 

equal  to  the  applied  strain,  £°.  This  is  also  equal  to  the  strain  without  cracks,  enc,  by  definition. 
The  average  strain  in  each  element  is  then: 


Sa-£  +£a  (3  8) 

where  £ad  is  the  disturbance  strain  in  each  element,  by  definition  the  difference  between  the  aver¬ 
age  overall  strain  and  the  average  strain  in  each  element.  The  disturbance  strains  in  the  a-element 
are  calculated  from  the  eigenstrains  in  the  a-element,  plus  a  volume  weighted  average  of  the  other 
eigenstrains.  In  this  case,  each  element  is  homogenized  by  the  overall  modulus;  that  is,  each  ele¬ 
ment  is  taken  to  have  the  overall  modulus  and  eigenstrains  are  introduced  to  keep  the  stresses  con- 
sistentUsing  that  assumption,  the  average  strains  and  stresses  are  found  for  each  element.  The 
modulus  is  defined  as  the  relation  between  the  average  stress  and  the  average  strain  for  the  entire 
RVE.  These  equations  are  solved  iteratively  to  reach  a  solution.  The  mathematical  details  are 
shown  in  Table  3.1. 

The  model  for  the  RVE  which  includes  cracks  follows  the  same  form.  The  cracks  in  this  case 
increase  the  strain  in  each  element,  but  make  no  contribution  the  stress.  The  additional  strain 
caused  by  the  cracks  is  found  by  calculating  the  crack  opening  displacement  (COD)  for  each 
crack.  The  average  stress  in  each  element  is  found  during  each  iteration  while  solving  the  prob¬ 
lem.  These  average  stresses  are  assumed  to  be  farfield  stresses,  and  the  COD  (and  therefore  the 
strain)  are  found  directly.  Using  a  method  similar  to  that  used  for  the  non-cracked  RVE,  the  aver- 
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age  stresses  and  strains  in  the  RVE  are  found,  and  the  overall  modulus  is  calculated.  The  equa¬ 
tions  are  solved  iteratively  to  find  the  solution. 

The  growth  of  the  cracks  is  based  upon  a  dynamic  crack-growth  model.  In  simple  terms,  a  crack 
tip  experience  a  crack  intensity  factor  Kj.  When  Kj  reaches  some  critical  crack  intensity,  Kj^,  the 
crack  will  extend.  So  it  is  necessary  to  find  both  Kr,  likely  using  linear  elastic  fracture  mechanics 
as  a  basis,  and  KIC,  which  will  take  into  account  the  fracture  properties  of  concrete  and  include  a 
dynamic  component.  Kj  is  a  function  of  the  material  properties  and  geometry,  as  well  as  the  gen¬ 
eral  loading  condition,  and  is  found  using  principles  of  linear  elastic  fracture  mechanics  (LEFM). 
While  the  limestone  and  mortar  inclusions  must  be  ellipsoidal,  the  cracks  can  be  any  geometry- 
penny  shaped,  ellipsoidal,  or  wing  cracks  are  all  possible. 

The  concrete  will  eventually  fail.  In  order  to  best  predict  the  failure,  the  size  of  the  cracks  are 
tracked  and  this  information  is  used  to  help  predict  the  failure.  The  sample  will  fail  when  KIC 
remains  less  than  Kr  By  tracking  the  crack  the  crack  length,  the  current  damaged  state  of  the 
material  can  be  assessed. 

Since  the  different  cracks  will  grow  at  different  rates,  an  anisotropy  may  develop  in  the  concrete. 
Even  before  any  loading  has  taken  place,  the  properties  in  compression  are  far  different  than  the 
properties  in  tension.  The  induced  anisotropy  and  the  difference  between  compressive  and  tensile 
response  are  both  addressed. 

In  order  for  this  work  to  be  useful,  it  must  be  able  to  be  discretized  for  use  in  numerical  simula¬ 
tions,  such  as  FEM  programs.  For  that  reason,  the  mathematical  complexity  in  the  final  model 
should  be  manageable.  Of  course,  the  most  important  feature  of  the  model  will  be  its  ability  to 
accurately  predict  the  response  and  failure  of  the  concrete.  The  model  is  a  physically  based,  rate 
dependent,  damage  model  that  can  be  used  in  simulations  to  accurately  predict  the  response  and 
failure  of  concrete.  The  calculation  of  modulus  for  the  cracked  and  non-cracked  RVE  have  been 
found  using  Mathematica,  and  the  code  is  included. 
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Table  3.1:  Mathematical  summary  of  model 


Quantity  |  Non-cracked  RVE 


Applied 

Strain 


Overall 

Modulus 


Nominal 

Stress 


Non-crack  nc  o 
s  =  s 

strain  _ 


7 77  NC  ,  d 

Strain  in  a  sa  =  e  +  £a 


Cracked  RVE 


0  ^NC  NC 

cr  =  C  :  £ 


0  _  77  O 

cr  =  C  :  s 


gNC  =  DNC  :  <x°  =  DNC  :  C":  e° 


sa  =  eNC  +  £a  +  £Q,crack 


Average 

Disturbance 


Definitions 


Crack  strain  |N/A 


Definitions 


eaA  =  Sa:  +  (Sv-  S„):  «/>*} 

a _ 


f  n 

€nd  =  SylSn*  +  /  — :  77  (Sv  -  Sa) :  {^or*  -  £P  | 

Z J  In  U  —  Id:}  B= 1 

_ q=l 

Sq,  :  EshelbyTensor  of  a  inclusion  (a  *  n) ,  Sy :  Eshelby  Tensor  of  outer  ellipsoid 

f  ;  Volume  fraction  of  a  £a  ■  Average  eigenstrain  in  a 

a  »  _ _ 


£*““*  =  Hq,'  :  (Ta=  Ha,':  Ca  :  (l(4SJ  +  Bg) :  :  q-u  =  Htt  : <ru 


Ha  =  Ha'  :Ca :  (l(4s)  +  Ba) :  DNC  Ba  :e  s  £a 

Ha’: Crack  relation  for  Dilute  distribution  model _ 


era  =  Ca  :  (sa-  S.™*)  =  Ca  :  (eNC  +  ej) 

=  Ctt : (1<4 s>  +  Ba) :  eNC  =  Ca :  (lc4 s)  +  Ba) :  DNC :  tr° 


~  /  cracK\  7=7  /  *  v 

Cqt  •  ^(x  J  —  C .  \&a  &a  ) 


Stress  in  a  cra-Ca:  sai=Ca  :(s  +ea 


Average 


Ca  :  (%)  =  C1^ :  (Sa  -  Sal _ 


Ca :  (fiNC  +  etf)  -  :  (eNC+  ead  -  €a* 


<  s  >=  eNC  +  <  sd  >=  (l*4s)  +  SV:A)  ic* 


<  £>=£NC+<£d>  +  <eaCraCk> 

=  (l(4s)+Sv:A+<H>  : C140)  : e1' 


A;  sht  =  <  i:*  >  :  Definition  of  concentration  tensor,  A .  <  ed  >=  Sv :  <£•*>,  Mori  -  Tanaka 


Definitions  Volume  average  <  (•)  >  =  7,  £»(•)« 

ct-l 


n 

<  <r  >=  2]f«  Ca  :  (%C  +£ad) 

0=1 

Average  =  («*+«.“- «.*) 

stress  »=> 

=  CNC :  (eNC  +  <  «d >  -  <«*  >) 

=  CNC :  (l(4s)  +  (Sy  —  l(4s>) :  A) :  £nc 


<  <r  >=  ^  fa  Ca  :  (£nc  +  £ad) 
a=l 

=  2jfa  CT:(eNC  +  Sad  -  Sa  +  £<7™^) 

0=1 

=  C:  («NC  +  <fid>  -  <£*>  +  <  CaCTaClt  >) 

=  C:(1(4s)  +  (Sv-1(4s)):A+<H>  :CNC):eNC 


C:  <  s>  =  <cr> 


Definition  of  ^nc  :<  e>  =  <cr> 

Modulus  _ _ _ 

C»==C>«:(1«  +  (Sv-1‘4->);A):{1«->+Sv:A)- 


Modulus 


Non-cracked: 
Cracked:  ^ 


C=C:(l(4s,  +  (Sv-l(4s,):A+<H>  :Cnc):(1(4s)+Sv:A+<H>  :Cnc) 
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4.0  Experiments 

4.1  Hopkinson  Bar  Tests  on  Concrete 

A  3-inch  Hopkinson  bar  is  used  to  test  the  concrete  at  high  strain  rates.  The  incident  and  transmit¬ 
ted  bars  are  each  8-feet  long.  All  the  bars  are  made  out  of  1040  steel,  and  elastic  wave  speed  is 
202,500  inches/  second.  This  steel  has  a  lower  hardness  than  the  maraging  steel  that  Hopkinson 
Bars  are  generally  made  from.  This  led  to  the  probability  of  the  bars  being  damaged.  To  prevent 
any  damage,  one-inch  thick  1040  steel  platens  are  placed  between  the  sample  and  the  bars. 

A  six  foot  brake  bar  was  used  to  stop  the  momentum  and  dissipate  the  energy  from  the  tests.  Two 
clamps  are  secured  around  the  brake  bar,  and  the  friction  between  the  clamp  and  the  brake  bar  is 
used  to  stop  its  momentum.  The  brake  is  what  stops  the  momentum  of  the  Hopkinson  Bar  and  dis¬ 
sipates  its  energy.  The  brake  works  by  clamping  a  three  inch  diameter  bar  (the  brake  bar)  that  is 
coaxial  with  and  next  to  the  transmitted  bar.  The  friction  between  the  clamp  and  the  brake  bar  dis¬ 
sipates  the  energy  in  the  system  as  the  brake  bar  slides  through  the  clamp.  A  torque  wrench  is  used 
to  tighten  the  clamp  the  desired  amount.  A  photograph  and  a  schematic  of  the  Hopkinson  bar  are 
shown  in  Figure  4.1. 


Brake,  Transmission  and 
Incident  Bars 


Figure  4.1:  Photograph  and  schematic  of  3-inch  Hopkinson  Bar 


Acquisition  equipment 

For  the  Hopkinson  Bar  tests,  similar  acquisition  equipment  was  used.  The  strain  gauges  were 
from  Micromeasurements  of  Raleigh,  North  Carolina,  and  were  attached  using  Loctite  Depend 
adhesive.  Two  gauges  are  located  on  every  incident  and  transmitted  bar,  and  are  used  in  a  Wheat¬ 
stone  Bridge  to  measure  the  strain,  using  30V  non-switching  DC  power  supplies  from  Kepco. 

The  digital  oscilloscope  used  to  capture  the  data  is  a  computer  controlled  Gagescope.  The  Gage- 
scope  is  a  pair  of  hardware  cards  that  are  installed  in  a  PC.  This  equipment  is  able  to  sample  either 
two  channels  of  data  at  up  to  100MHz  or  4  channels  at  up  to  50MHz,  using  a  12-bit  digitizer.For 
the  present  work,  samples  rates  around  2MHz  are  used.  The  Gagescope  is  controlled  by  software 
that  runs  under  MS-DOS. 
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Loading  Profiles 

For  the  Hopkinson  Bar  experiments,  the  elastic  wave  profile  affects  the  material  response.  In  order 
to  maintain  a  constant  strain  rate  during  an  experiment,  a  squarish  wave  profile  is  needed.  If  the 
rise  time  is  too  quick,  the  sample  may  not  come  into  equilibrium,  and  the  test  results  may  be  sus¬ 


pect.  An  elastic  wave  with  a  slower  rise  time  is  more  likely  to  result  in  establishing  equilibrium 
during  a  test,  but  a  constant  strain  rate  cannot  be  maintained.  Figure  4.2  shows  two  examples  of 
what  the  incident  pulse  and  strain  rates  typically  were  during  experiments.  In  the  rest  of  this 


paper,  these  two  classes  will  be  broadly  referred  to  as  square  and  triangular  pulses. 


Figure  4.2:  Typical  elastic  wave  (incident  pulse)  profiles  and  strain  rates  used  in  testing 


The  square  waves  are  normally  generated  during  the  striker  bar-  incident  bar  impact.  In  order  to 
create  a  triangular  pulse,  a  damping  material  must  be  used.  For  these  experiments,  commercially 
pure  aluminum  was  used  to  aid  in  the  generation  of  the  triangular  pulse.  This  material  was  chosen 
due  to  its  uniformity,  and  its  low  rate  of  strain-hardening.  It  was  important  to  choose  a  material 
which  had  little  strain  hardening;  otherwise  the  incident  and  striker  bars  could  be  damaged. 

Uniaxial  tests 

Typical  unixial  samples  are  shown  in  Figure  4.3.  The  samples  were  all  3”  in  diameter,  and  ranged 
in  length  from  0.75”  to  2.25”.  During  the  tests,  the  samples  were  held  in  place  between  two  1” 
thick  platens  with  tape.  The  platens  were  secured  using  the  rings  seen  below. 
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Figure  4.3:  Typical  samples  used  in  unconfined  concrete  tests. 


Tests  were  conducted  using  both  triangular  and  square  wave  profiles  at  a  number  of  strain  rates. 
The  strain  rate  effect  can  be  seen  in  Figure  4.4.  The  samples  in  this  tests  were  each  1”  long,  and 
were  destroyed  after  the  test  by  the  reloading  of  the  sample.  The  strain  rates  reported  below 
roughly  correspond  to  the  peak  strain  rates,  and  the  results  for  the  triangular  and  square  wave  pro¬ 
files  cannot  be  directly  compared.  However,  it  is  observed  that  the  samples  with  the  triangular 
wave  profile  have  a  smoother  measured  response. 


Figure  4.4:  Strain  rate  effect  on  unconfined  concrete  using  triangular  (left)  and  square  loading  profiles 


Uniaxial  tests  with  high  speed  photography 

High  speed  photographs  were  taken  during  some  experiments.  Some  samples  were  photographed 
without  cutting  parallels,  as  shown  in  Figure  4.3.  Other  samples  had  parallel  surfaces  cut  across 
the  diameter  of  the  sample,  similar  to  what  is  shown  in  Figure  4.12.  Photographs  of  the  tests  are 
shown  in  Figure  4.5  to  Figure  4.8,  where  the  loading  is  horizontal.The  strain  rate  at  failure  was 
about  250/s.  Half  of  the  results  presented  show  an  entire  face,  and  the  others  focus  on  a  smaller 
section  of  the  sample.  A  Hadland  Imacon  792  was  used.  This  camera  uses  one  3”  x  5”  piece  of 
Polaroid  film,  and  8  images  were  taken.  These  images  can  be  taken  at  rates  from  of  50  thousand  to 
20  million  frames  per  second.  The  resolution  is  10  lines/mm,  which  results  in  an  image  size  of 
about  200  x  200  pixels.  For  these  tests,  the  frame  rate  used  was  50,000/s  (20ps). 
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Figure  4.6:  High  speed  photographs  of  uncut  section  of  concrete 
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Cracks  can  be  seen  growing  axially  and  radially  in  Figure  4.5.  In  this  test,  one  crack  became  dom¬ 
inant  as  it  grew  the  quickest  across  the  sample.  The  large  crack  grew  across  the  mortar,  but  other 
cracks  are  observed  which  split  some  limestone  aggregates. 

Figure  4.6  shows  lots  of  crack  formation  after  the  peak  stress  has  been  reached.  Some  cracks  at 
this  stage  are  growing  parallel  with  the  loading  axis,  but  the  most  active  cracks  are  perpendicular 
to  the  initial  loading  direction 

In  Figure  4.7  a  crack  can  clearly  be  seen  propagating  at  the  sharp  tip  of  a  limestone  aggregate. 
Under  careful  observation,  a  crack  can  be  seen  growing  from  the  tip  of  a  limestone  piece  in  Figure 
4.8.  These  pictures  illustrate  a  crack  growing  across  the  sample  from  right  to  left.  In  all  cases,  sig¬ 
nificantly  more  crack  growth  was  observed  after  the  failure  stress.  However,  the  resolution  of  the 
camera,  as  well  as  the  total  number  of  frames,  limited  what  could  be  observed. 

Triaxial  tests 

Stainless  steel  T-bolt  band  clamps  (similar  in  shape  to,  but  stronger  than,  hose  clamps)  provided 
the  radial  confinement  for  some  triaxial  tests.  The  clamps  have  a  thickness  of  0.025”,  and  are  each 
0.75"  wide.  The  clamps  are  wrapped  around  the  samples  and  loaded  to  a  specified  confining  pres¬ 
sure.  A  thin  sheet  of  Teflon  is  located  between  the  clamp  and  the  sample  to  reduce  friction  in  the 
axial  direction.  On  some  samples,  platens  were  clamped  along  with  the  concrete;  on  the  other 
samples  the  platens  were  taped  to  the  clamps.  In  all  cases,  the  plastic  tubes  held  the  platens  (and 
samples)  in  place.  A  torque  wrench  was  used  to  load  the  samples. 


Figure  4.9:  Topical  confined-clamp  sample 


On  a  number  of  samples,  the  relation  between  the  applied  torque  and  the  radial  confinement  was 
calibrated;  the  hoop  strain  in  the  clamp  was  measured  using  a  strain  gage.  The  radial  stress  and  is 
calculated  using  that  measurement.  On  one  sample,  the  hoop  strain  in  the  clamp  was  measured 
dynamically.  In  this  case,  the  radial  and  volumetric  strain  in  the  sample  were  calculate. 

As  shown  in  Figure  4.10,  the  steel  band  clamp  clearly  yielded  possibly  before  and  certainly  dur¬ 
ing  the  test.  This  suggests  that  after  an  initial  rise  in  the  confining  pressure,  it  increases  much 
more  slowly.  Thus  an  initial  difference  in  confining  pressure  is  unlikely  to  produce  different 
results. This  test  was  done  at  a  nominal  strain  rate  of  200/s  with  a  triangular  wave. 
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Figure  4.10:  Dynamic  measurement  of  hoop  and  volumetric  strain.  The  figure  on  the  left  shows  that  the  steel 
clamp  yielded  during  the  experiment;  the  figure  on  the  right  calculates  the  concrete  response  during  the  test 


As  shown  in  Figure  4.1 1,  the  two  samples  were  initially  subject  to  different  confining  pressures. 
However,  the  response  of  the  two  samples  are  nearly  identical,  due  to  the  increase  in  confinement 
during  the  test.  The  samples,  each  1.5”  long,  were  tested  at  a  strain  rate  of  about  200/s,  with  a  tri¬ 
angular  elastic  wave  profile. 


Figure  4.11:  Effect  of  initial  confinement  on  concrete  response  and  response  of  radially  confined  and  uncon- 
fined  concrete. 

The  presence  or  absence  of  confinement  has  a  dramatic  effect  on  the  material  response  and  fail¬ 
ure.  Figure  4.11  shows  an  example  for  a  pair  of  2  inch  long  specimens  which  have  been  tested  at  a 
strain  rate  of  200/s  using  a  square  loading  profile.  After  testing,  confined  samples  appeared  largely 
undamaged  while  the  unconfined  samples  have  been  shattered.  One  reason  the  unconfined  sam¬ 
ples  were  turned  into  rubble  is  because  they  were  subject  to  repeated  loading  by  the  Hopkinson 
Bar.  Micrographs  taken  of  the  confined  samples  after  testing  are  presented  in  the  next  section. 


Microscopy  on  damaged  samples 

Most  of  the  samples  tested  with  the  clamp  confinement  remained  largely  intact  after  the  experi¬ 
ments.  Some  of  these  samples  were  used  for  microscopy,  and  the  micrographs  are  presented  here. 
The  cross  section,  one  inch  thick,  that  was  examined  is  shown  in  Figure  4. 12.  Both  top  and  bottom 
were  cut  with  a  diamond  saw,  and  one  face  was  polished. 
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Concrete  Sample 


Figure  4.12:  Cross  section  used  for  microscopy 

The  first  pictures  were  taken  of  a  sample  that  had  been  loaded  with  a  triangular  loading  profile  at 
a  peak  strain  rate  of  200/s,  under  200psi  initial  confinement,  sample  cnc3-2q-4.  The  arrows  show 
the  direction  of  the  axial  straining.  An  axial  crack  through  a  rock  is  plainly  visible  in  Figure  4.13. 
Sample  cnc3-2q-2  was  2.25"  long  and  loaded  with  a  trapezoidal  pulse  at  a  strain  rate  of  100/s  with 
an  initial  confinement  of  lOOOpsi.  The  sample  was  tested  a  total  of  3  times.  The  two  photographs, 
Figure  4. 14,  focus  on  what  appears  to  be  a  crack  running  part  of  the  way  through  a  rock  in  the 
sample.  Most  of  the  rocks  did  not  have  visible  cracks  in  the  axial  direction. 


Figure  4.13:  Axial  crack  in  sample  cnc3-2q-4,  at  5x  (left)  and  20x  (right) 
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Figure  4.14:  Crack  in  cnc3-2q-2, 5x  (left)  and  20x  (right) 


Pneumatic  Pressure  Vessel 

A  large  diameter  pressure  vessel  provides  the  pneumatic  confinement;  see  Figures  3  and  4.  A  steel 
cylinder  of  6  inches  outside  diameter  and  half-inch  wall  thickness  surrounds  the  sample.  The  alu¬ 
minum  end  plates  are  each  1  inch  thick,  and  are  held  together  by  four  3/8  inch  bolts.  A  rubber 
sleeve,  secured  by  hose  clamps,  is  placed  over  the  sample  to  prevent  direct  contact  between  the 
sample  are  the  high-pressure  gas.  O-rings  are  used  to  seal  the  pressurized  chamber.  A  band  clamp 
and  the  brake  bar  restrain  the  incident  and  transmission  bars,  respectively,  to  prevent  them  from 
moving  apart.  This  pressure  vessel  provides  a  constant  radial  confinement  of  up  to  1000  psi. 
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4.2  Experiments  on  rock  and  mortar 

In  addition  to  testing  the  concrete,  the  constituent  materials  have  also  been  tested.  Ultrasonics 
have  been  used  for  NDE,  an  Instron  machine  has  been  used  for  low  strain  rate  uniaxial  testing, 
and  Hopkinson  bars  have  been  used  for  uniaxial  and  triaxial  tests  at  high  strain  rates. 


Low  strain  rate  tests  on  mortar 

Mortar  samples  prepared  without  the  limestone  aggregate  have  been  tested.  Three  samples,  each 
about  0.75”  long  and  0.7”  in  diameter,  were  tested  at  low  strain  rates.  A  20kips  hydraulic  Instron 
machine  is  used  for  mnning  low  strain  rate  experiments.  National  Instruments  hardware  and  soft¬ 
ware  (Labview)  were  used  for  the  data  acquisition.  Data  was  sampled  at  rates  from  1/s  to  10/s.  In 
the  tests,  the  load  measurement  is  very  accurate.  However,  the  strain  measurement  is  not  as  accu¬ 
rate.  In  both  cases,  the  strain  measurement  is  largely  dependent  on  how  well  the  strain  gage  or 
extensometer  is  attached  to  the  sample.  The  attachment  of  the  strain  gage  to  the  mortar  is  good  at 
low  strains,  but  at  higher  strains  the  deformation  may  be  affected  locally  by  the  strain  gage  and  the 

glue  used  to  attach  it.  At  low  strains,  Young’s  Modulus  is  found  to  be  4.7xl06  psi.  One  result  is 
shown  in  Figure  4.16. 


High  strain  rate  uniaxial  tests  on  mortar 

Uniaxial  tests  were  done  on  the  mortar  samples  using  the  3/4”  recovery  Hopkinson  bar.  The  3/4” 
Hopkinson  bar  is  used  to  test  the  rock  and  mortar  at  high  strain  rates.  This  Hopkinson  Bar  is  made 
from  Maraging  steel,  and  the  elastic  wave  velocity  was  measured  as  200,000  in/s.  The  incident 
and  transmitted  bars  are  each  48”  long,  and  striker  bars  between  3  and  18  inches  may  be  used  with 
this  tool.  The  bar  is  designed  to  prevent  reloading  of  the  sample,  and  the  use  of  the  recovery  Hop¬ 
kinson  bar  allows  for  valid  post-experiment  assessment  of  samples.  This  assessment  includes 
NDE  with  ultrasonics,  or  ultimate  strength  tests  using  low  strain  rate  testing  equipment.  As  shown 
in  Figure  4.16,  the  small  difference  in  strain  rates  had  little  effect  on  the  response,  but  more  tests 
are  needed  to  verify  this.  On  the  right,  the  response  of  the  mortar  at  high  and  low  strain  rates  is 
contrasted.  It  is  obvious  that  the  failure  stress  is  much  higher  at  the  high  strain  rates. 
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Strain  (%) 


Figure  4.16:  Rate  effects  on  mortar  response 
High  strain  rate  triaxial  tests  on  mortar 

The  classical  Hopkinson  Bar  can  be  modified  to  allow  dynamic  triaxial  compressive  loading  of  a 
sample.  This  technique  simultaneously  loads  the  sample  in  the  axial  and  radial  directions.  Figure 
4.17  and  Figure  4.18  are  a  photograph  and  schematic  of  this  system.  The  striker  (A)  impacts  the 
first  incident  bar  (B),  generating  the  incident  pulse.  The  wave  is  transmitted  to  the  second  incident 
bar  (C)  and  the  incident  tube  (F).  (Note:  With  this  design,  the  reflected  pulse  can  not  be  measured 
directly,  so  it  is  calculated  as  the  difference  between  the  transmitted  pulse  and  the  incident  pulse.) 
The  sample  (D)  is  inside  a  Teflon  tube  (G),  which  in  turn  is  inside  an  aluminum  sleeve  (H). 


Figure  4.17:  Photograph  of  triaxial  load  cell 

The  confinement  is  provided  by  the  Teflon,  which  is  dynamically  compressed  between  the  inci¬ 
dent  tube  (F)  and  the  transmission  tube  (I).  Restrained  laterally  by  the  aluminum  sleeve,  a  large 
hydrostatic  stress  is  produced  in  the  Teflon.  This  pressure  creates  a  large  radial  stress  on  the  sam¬ 
ple.  The  hoop  strain  in  the  aluminum  sleeve  is  measured,  and  the  radial  confining  stress  is  calcu¬ 
lated.  The  radial  stress  can  be  controlled  independently  from  the  axial  stress  and  strain  to  a  limited 
extent,  by,  e.g.,  altering  the  thickness  of  the  aluminum  sleeve  to  control  when  the  sleeve  yields. 
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The  simultaneous  loading  in  the  radial  and  axial  directions  is  assured  by  the  design  of  the  bar-  the 
stress  waves  in  the  incident  bar  and  incident  tube  are  generated  at  the  same  time,  they  are  made  of 
the  same  material,  and  they  have  nearly  the  same  length.  Thus,  the  stress  wave  in  the  incident  bar 
reaches  the  sample  (loading  it  axially)  at  the  same  time  as  the  stress  wave  in  the  incident  tube 
reaches  the  Teflon  (loading  the  sample  radially). 


This  method  has  been  used  to  test  several  samples.  The  confinement  had  a  tremendous  effect  on 
the  response  of  the  mortar.  Figure  4.19  shows  the  response  during  three  tests.  Samples  TXMK  and 
TXML  were  tested  under  identical  conditions,  except  TXML  used  the  thinner  aluminum  sleeve. 
During  that  test,  the  sleeve  yielded,  which  resulted  in  a  lower  confining  pressure,  which  is  esti¬ 
mated  in  that  figure.  The  lower  radial  stress  led  to  lower  axial  stresses.  Comparison  with  Figure 
4.16  shows  the  dramatic  effect  confinement  has  on  the  mortar  response. 


Figure  4.19:  Response  of  mortar  under  triaxial  compression,  triaxiai  pieces  un-assembied. 

A  scanning  electron  microscope  (SEM)  has  been  used  to  observe  some  of  the  tested  an  untested 
samples.  Figure  4.20  shows  a  micrograph  of  an  untested  concrete  sample  and  a  tested  mortar  sam- 
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pie.  Many  more  cracks  are  observed  in  the  tested  samples  than  the  untested  samples,  when 
observed  under  similar  magnifications. 


Figure  4.20:  SEM  micrographs  of:  untested  concrete  (25x);  mortar  tested  in  triaxial  compression  (500x) 


1  ypw  strain  rate  tests  on  limestone 

Low  strain  rate  tests  were  done  using  the  20  kips  Instron  machine.  The  samples  were  0.5”  in 
diameter  and  0.5”  long.  The  data  shows  that  the  rocks  fail  elastically  under  uniaxial  compression. 
As  is  the  case  with  the  mortar,  the  load  measurement  is  very  accurate.  The  displacement  measure¬ 
ment  should  be  more  accurate  than  it  is  for  the  mortar,  because  the  strain  gages  adhere  much  bet¬ 
ter  to  the  surface  of  the  limestone  than  they  do  to  the  mortar.The  measured  modulus  for  the  rock  is 

about  10.5  x  106  psi.  See  Figure  4.21. 

High  strain  rate  uniaxial  tests  on  limestone 

The  3/4”  Hopkinson  bar,  was  used  to  conduct  uniaxial  tests  on  the  limestone.  The  use  of  the 
recovery  Hopkinson  Bar  allowed  for  the  retrieval  of  damaged  by  intact  samples,  which  are  evalu¬ 
ated  using  NDE. 


Figure  4.21:  Comparison  of  rock  and  mortar  response  (left);  high  and  low  strain  rate  limestone  response 
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Figure  4.21  shows  the  response  of  the  limestone.  Comparison  with  the  results  from  low  strain  rate 
testing  show  that  limestone  has  a  greater  failure  stress  under  dynamic  conditions.  However,  due  to 
the  variability  from  sample  to  sample,  no  definite  conclusion  can  be  drawn  at  this  time,  though  it 
is  suspected  that  this  is  a  real  effect.  Above,  Figure  4.22  shows  a  comparison  of  the  rock  and  mor¬ 
tar  at  high  strain  rates.  The  limestone  clearly  has  a  greater  failure  stress  than  the  mortar. 


Figure  4.22:  Comparison  of  rock  and  mortar  response  at  high  strain  rates 
NDE  of  limestone  and  mortar  using  ultrasonics 

Ultrasonic  equipment  can  be  used  for  NDE  of  a  range  of  materials.  Samples  may  be  evaluated 
with  ultrasonics,  tested  on  a  recovery  Hopkinson  Bar,  then  tested  again  with  ultrasonics  to  assess 
changes.  The  equipment  used  in  the  current  study  is  described  briefly  below. 

The  wave  generation  and  initial  transmission  to  the  sample  is  identical  to  the  procedure  used  by 
Aashish  Rohatgi  in  his  research  at  UCSD.  He  describes  the  set-up,  in  part,  as  follows  [Rohatgi, 
pp.42-45]: 

A  constant  sinusoidal  wave  (CW)  of  the  desired  frequency  and  amplitude  (~125mv)  is  produced  by  the  Marconi 
Synthesizer  and  supplied  to  the  waveform  generator  in  the  MBS-8000  ultrasonic  measurement  system,  which  in 
turn  generates  a  high  voltage  (-100  V)  signal  (RF)  of  the  same  frequency.  The  pulse  width,  amplitude  and  fre¬ 
quency  of  the  RF  signal  are  controlled  via  a  computer  interface.  The  RF  pulse  excites  the  broad  band  transducer 
which  sends  a  high  frequency  pulse... 

The  transducer  is  put  into  contact  with  the  sample,  using  water  as  a  couplant  for  longitudinal  tests, 
and  honey  as  a  couplant  for  shear  test.  Another  transducer  is  placed  on  the  opposite  end  of  the 
sample,  which  is  excited  by  the  acoustic  wave  and  generates  a  voltage,  which  is  in  then  measured. 
Comparison  of  the  initial  signal  and  the  transmitted  signal  allows  for  computation  for  the  velocity 
through  the  sample. 

The  acquisition  equipment  is  the  same  digital  Gagescope  that  is  described  in  section .  In  this  case, 
the  Gagescope  samples  the  data  at  60MHz. 

Transducers  with  natural  resonant  frequencies  of  5MHz  and  10MHz  were  used.  The  transducers 
were  supplied  by  Panametrics.  Both  shear  and  longitudinal  transducers  were  used  in  the  experi¬ 
ments,  and  all  contained  a  piezoelectric  crystal. 
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Measurements  of  the  acoustic  wave  velocities  were  done  on  both  the  rock  and  mortar  samples, 
plus  maraging  steel  for  comparison.  The  elastic  modulus  and  the  Poisson  ratio  can  be  calculated 
from  the  plane  (Ci)  and  shear  (C2)  wave  velocities,  when  the  density  is  also  known. 


_  Cf  -  2  x  Cl 

V“2  X(C?-C?)  (4.1) 

C12xpx(l  +  v)x(l-2v) 

('-V)  (4.2) 

Inspection  of  the  calculated  elastic  constants  for  steel  reveal  that  the  method  is  valid.  The  differ¬ 
ence  between  the  elastic  modulus  found  via  low  strain  rate  testing  (4.7  x  106  psi)  differs  signifi¬ 
cantly  from  the  ultrasonic  testing  (3.3  x  106  psi).  The  low  strain  rate  testing  will  be  repeated, 
taking  care  in  the  displacement  measurement,  which  can  be  difficult  to  make  accurately. 


Table  4.1:  Wave  velocities  and  calculated  elastic  constants  of  rock  and  mortar 


Material 

Specific 

Density 

Wave  Velocities  (103  in/s) 

Poisson 

Ratio 

Young’s  Modulus 
(106  psi) 

Range 

Plane  (Cj) 

Shear  (C2) 

Rock 

2.69 

258 

134 

0.32 

11.9 

2.5% 

Mortar 

2.11 

131 

85 

0.15 

3.3 

1.6% 

Maraging  Steel 

8.06 

223 

117 

0.31 

27.2 

N/A 

The  results  from  the  ultrasonic  testing  are  summarized  in  Table  4.1.  One  test  result  is  shown  in 
Section  4.23.  For  the  limestone,  the  two  separate  calculations  of  the  elastic  modulus  are  quite 

close,  10.5  x  106  psi  for  the  low  strain  rate  testing,  and  1 1.9  x  106  psi  for  the  ultrasonic  testing. 
Some  of  these  tests  may  be  repeated  to  ensure  a  statistically  sound  database. 
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Figure  4.23:  Measurement  of  the  shear  wave  speed  in  rock  using  ultrasonics 


5.0  Material  characterization 

Examples  have  already  been  shown  of  optical  and  SEM  micrographs  used  to  aid  in  the  character¬ 
ization  of  the  sample.  Those  pictures  yield  sound  information  on  die  size  and  density  of  cracks  in 
the  materials.  To  measure  the  volume  fraction  of  limestone,  a  direct  quantitative  method  is  used. 
First,  a  concrete  sample  is  scanned  using  a  computer  scanner,  Figure  5.1.  Then,  the  limestone 
aggregates  are  outlined  with  the  aid  of  a  graphics  tablet.  The  limestone  is  then  set  to  black,  and  the 
software  counts  the  percent  of  black  pixels.  In  the  example  shown  below,  the  limestone  has  a  41% 
volume  fraction.  This  information  is  then  used  to  determine  the  size  of  the  inclusions  in  the 
model. 


Figure  5.1:  Scan  of  concrete  sample  with  limestone  outlined;  scan  of  concrete  with  limestone  darkened. 


6.0  Conclusions 

A  RVE  has  been  designed  to  embody  the  primary  features  of  the  concrete,  based  upon  its  physical 
microstructure.  This  model  is  designed  using  information  from  material  characterization  and 
experimental  results.  Using  micromechanics,  fracture  mechanics  and  damage  mechanics,  the 
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response  and  failure  of  the  concrete  is  predicted.  The  modulus  is  found  by  numerically  imple¬ 
menting  the  model  in  Mathematics  The  crack  growth  and  sample  failure  criteria  are  being  added 

to  the  current  implementation. 
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ADDENDUM  TO  FINAL  REPORT  -  Mathematica  file 


double _crack_6.nb 


(*  Set-up  problem  *) 

(*  Physical  parameters  *) 

Clear [ "Global' * " ] ;  (*  clears  all  variables  from  memory  *) 
n  =  2 ;  (*  Number  of  micro-elements  *) 

Vref  =  1/3;  jUref  =  1;  (*  Initial  properties  of  body  *) 

Kref  =  2  Uref  (1  +  Vre£)  ;  (*  Initial  Bulk  Modulus  of  Body  *) 

3(1-2  vref ) 

/il  -  50;  vx  =  1  /  3 ;  ju2  =  1;  v2  =  1  /  3;  (*  Properties  of  micro-elements  *) 


Youngs 2  =2  \x2  (l  +  v2);  (*  Young'  s  modulus  of  mortar  *) 


fcrackx  =  0?  fcrack2  =  0;  fcrack3  =  .1; 


fi  =  .4;  f2  =  1  -  f i ;  (*  Volume  percent  of  micro-elements  *) 

nclterations  =  10; 

Iterations  =  25;  (*  Number  of  loops  to  iterate  over  *} 


weighty  =  0; 


weightn 


(1  -  (1  -  f„)  weight DOtn) 


Print  [  {weightn  ,  weighty }  ]  ? 

(*  Define  basic  functions  used 
to  create  and  transform  elasticity  and  Eshelby  Tensors  *) 

(  a  The 

matrix  needed  to  convert  the  compliance  tensor  to  a  matrix,,  and  its  inverse. 


W  =  Table  [  0 ,  {il,  6},  {i2,  6}]; 


Do [w[ [i, 

i]]  =  1. 

<i.  3}]; 

Do [ W [ [ i , 

i]  ]  =2, 

{i.  4,  6}] 
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InvW  =  Inverse  [W]  ; 


Sec.  up 


Lze  and  mid  ax  values 


ior  some  vectors  and  matrices  *} 


Zero  6  =  Table  [  0 ,  {i,  6}];  {*  Length  6  vector  with  a  magnitudes  *) 


Zero 6by 6  =  Table  [  0,  {il,  6},  {i2,  6}]; 

A  =  Zero6by6;  Aprime  =  Zero6by6;  Aprimenew  =  Zero6by6; 

Acrack  = Zero6by6; 

distrainl = Zero6by6; 

distrain2  = Zero6by6; 

Do  [  eigeng  =  Table  [xSi<  {i,  6}],  {0,  n}  ]  ;  (*  Fills  out  eigenstrain  vectors  *) 

{*  widgets  to  populate  and  transform  elasticity  and  eshelby  tensors  *) 
first [1]  = 1;  first [2]  =  2;  first [3]  =  3;  first [4]  =  2;  first [5]  =3;  first [6]  = 1; 
second [  1  ]  =  1;  second[2]  =2;  second[3]  =3;  second[4]  =3;  second[5]  =1;  second[6]  =2; 


MakeTwo  [Mod4_]  :  = 

Module  [{DumMod2},  DumMod2  =  Table  [  0 ,  {il,  6),  {i2,  6}  ] ;  Do[DumMod2  [  [i,  j  ]  ]  = 

Mod4  [  [first  [i]  ,  second[i],  first[j],  second[j]]],  {i,  6},  {j,  6}  ];  Return  [DumMod2]  ;  ] 


PopEshSphe  re [ vDum_]  :  = 

Module  [  {DumMod} ,  DumMod  =  Table  [ 0 ,  {il,  3},  {i2,  3},  {i3,  3},  {i4,  3}]; 


Do  [ DumMod [  [i,  j,  i,  j]] 
Do  [DumMod  [  [i,  j,  j,  i]  ] 
Do  [DumMod  [  [i,  i,  j,  j]] 
Do  [DumMod [  [i,  i,  i,  i]] 


(4-5  vDum ) 

15  (1  -  vDum)  ' 
(4-5  vDum) 

15  (1  -  vDum)  ' 
5  vDum  -  1 
15(1-  vDum)  ' 
/  5  vDum  -  1 

\  15  (1  -  vDum) 


{i,  3},  {j,  3}]; 
{i,  3},  {j,  3}]  ; 

{i,  3},  {j,  3}]  ; 

2  .(4  -  5  vDum)  \ 

+  15  (1  -  vDum)  j  ' 


{i.  3}] 


Return [DumMod] ; ] ; 


PopIsoMod[juDum__,  vDum_]  :  = 

Module  [  {DumMod}  ,  DumMod = Table [ 0 ,  {il,  3},  {i2,  3},  {i3,  3},  {i4,  3}]; 
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{*  Initial  guess  for  shear  modulus  of  matrix  *) 


Cref  =  MakeTwo  [ PopIsoMod  [p/M /  vref  ]  ] ;  {*  Initial  Reference  elasticity  *} 
(*  Check  for  relevance  of  crack  strain  *} 

(*  Define  functions  specific  for  this  problem  *) 

(*  Eigenstrains  averaged  over  non-or  elements  *) 


0 1 he rEigen Average  [a_,  n2_]  :=  Module  [  {Result} ,  Result  -  Table  [0,  {il,  6},  {i2r  6}]? 

n2  j; 

Result  =  (  1  -  eigen^  J  -  ±  _g£  eigena;  Return  [Result  ] ;  ]  ; 

0=1  a  01 

(******************************************************************************** 

Do[{ 

(*  Define  and  evaluate  Disturbance  strain  for  all  micro-elements  *) 

StrainDisturb  [a_,  nl_]  :  =  Module  [  (Result } ,  Result  = Table [ 0 ,  {il,  6},  {i2,  6}]; 

Result  =  8a  .  W  *  eigena  +  (sy  -  8a)  . W  .  Other EigenAverag e  [a,  nl] ;  Return  [Result  ];  ] 
(*  defines  disturbance  strain  for  all  microelements  except  outer  ellipsoid  *) 

Do [edista  |  StrainDisturb [a,  n]  ,  {a,  n-1}];  (*  Evaluate  disturbance  strain  *) 

Edistn  =  Sv  .  W  .  eigenn  +  \  j-  g-y  (S*  -  Sa )  .  W  .  eigena  -  £  f,  eigen0  ; 

a=l  ^ 

(*  Find  disturbance  strain  for  outer  ellipsoid  *) 

Do[{ 

(*  Sets  up  each  case  for  one  non-zero  component  *) 

Zero6;  (*  resets  e  to  zero  *) 
cigenew1  = Zero6;  cigenew2  =  Zero 6; 
eigenAve  =  Zero 6; 

e[[y]]  =  l/lntegerPart[l  +  yT-];  (*  Set  one  component  of  e=l  '*) 

(*  Solve  for  eigenstrain  for  each  case  *) 

{ { { e igenew1 ,  eigenew2}}}  = ReplaceAll [ {eige^ ,  eigen2}, 

{Solve  [  {Ci  .  W  .  (e  +  edisti)  -  Cre£  •  W  .  (e  +  edisti  -  eigenj )  ==Zero6, 

C3  .W.  (e  +  Edist2 )  -Cref  .W.  (e  +  edist2  -sigen2)  ==Zero6>, 

{ (Xl }  2  ,  (Xi)2.  (Xi  )  3  /  (Xi)4, 

{Xi  )  5  ,  (Xi  )  6  ,  (X2  )  !  /  (X2  )  2  /  (X2  )  3  /  {X2  )  4  /  {X2  }  5  /  (X2  )  g  }  ]  }  ]  ' 
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(*  Finds  volume-average  eigensirain  *) 

n 

eigenAve  =  eigenewa; 

a=l 

Do [A [ [ i ,  y]  ]  =  eigenAve [ [i]  ] ,  {i,  6}] ; 

},  {r.  6}]; 

C„e„  =  Cret  .  W  .  (InvW  +  (Sv  -  InvW)  .  W  .  A)  .  Inverse [ InvW +  Sv  .  W  .  A]  .  InvW; 

Cref  =  ^new  r 

/^new  =  Cref  [  [  4  ,  4  ]  ]  i 

,  Cref  [[1,  2]] 

VneW  '  2  (Cref  [[1,  2]]  +  f-J-new )  / 

„  2  (J-new  (  1  +  ^new  ) 

Kne„-  3{1_2vnew)  ' 

Clear [StrainDisturb] ; 

Sv  = MakeTwo [ PopEshSphere [ vnew ] ] ;  {*  Eshelby  tensor  of  outer  ellipsoid  *} 

Do [sa =  MakeTwo [ PopEshSphere [vn ] ] ,  {a,  n-1}];  (*  Eshelby  tensor  of  micro-elements  *) 


]  ,  {nc Iterations} ] ; 


{Dcrack2,  C crack  }  =  PennyCrackSt ress  [ Cref  ,  vn,  jun  ,  fcracki ,  fcrack2 ,  f  crack3  ]  ; 


C  nocrack  —  Cref  1 

Dnocrack  =  IHVW  .  Inverse  [Cnocraok  ]  .  InvW; 

(*  Find  the  concentration  tensor  *) 

(*  Large  do-loop  *) 

Do[{ 

(*  Define  and  evaluate  Disturbance  strain  for  ail  micro-elements  *) 

StrainDisturb [a_,  nl_J  := Module [ {Result} ,  Result  =  Table  [ 0 ,  (il,  6},  {i2,  6}]; 

Result  =  sa  .  W  •  £igena  +  (Sv  -  Sa)  .  W  .  OtherEigenAverage  [a,  nl] ;  Return  [Result  ] ;  ]  ; 
(*  defines  disturbance  strain  for  all  microelements  except  outer  ellipsoid  *} 
Do[edista  =  StrainDisturb [a,  n]  ,  {a,  n-1}];  (*  Evaluate  disturbance  strain  *) 

l  n 

edistn  =  Sv  .  W  .  eigenn  +  \  ^  j-y  (Sv-Sa)  .  W  .  eigen^  -  ^  f e  eigen^  ; 

0=1  '  ^-1 

{*  Find  disturbance  strain  for  outer  ellipsoid  *) 

Do[{ 

(*  Sets  up  each  case  for  one  non-zero  component  *) 

e  =  Zero6;  (*  resets  e  to  zero  *) 

eigenew1  =  Zero 6;  eigenew2  =  Zero6;  edistnewi  =  Zero6;  sdistnew2  =  Zero6; 
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sigenAve  =  Zero  e'¬ 
er  [[y]]  =  l/integerPart  [l  +  j^-]  ;  {*  Set  one  component  of  e=l  *) 

Dcrack  -  Dcrack2  .  W  .  C2  .  W  .  (InvW  +  Aprimenew)  .  W  •  Dno  crack  ? 
ecrackref  —  Dcrack  •  W  •  CJlocracjc  »  W  •  £ ; 
scrackn  =  weightn  x  ecrackref ; 

Do[ecracka  =  weightnotn  x  ecrackref ,  { a ,  n-1}]; 


(*  Solve  for  eigenstrain  for  each  case  *) 

{{{eigenew1(  eigenew2,  Edistnewj ,  edistnew2 } } }  = 

ReplaceAll[{eigen1(  eigen2,  edisti,  edist2}, 

{Solve  [{Ci  .W.  ( e  +  edistj )  -  CreE  .  W  .  (e  +  ccracki  +  cdistj  -  eige^ )  ==  Zero6, 

C3  .  W  .  (e  +  edist2 )  -Cre£  .  W  .  (e  +  ecrack2  +  Edist2  -  Eigen2)  ==  Zero6}, 

{(Xl)ir  (Xl)2,  (Xi)3, 

(Xi  )  4  /  {Xi  )  j  ,  {Xj  )  4  ,  (X2  )  1  ,  (X2  )  2  t  (X2  )  3  t  {X2  )  4  t  (X2  j  5  r  {x2  )  4  }  ]  }  ]  t 

(*  Finds  volume -average  eigenstrain  *) 

edistfinali  =  Zero6; 

edistfinal2  =  Zero6; 

(*  Junk=ReplaceAll [edisti , 

{(Xi)^,  (Xi)2,  (Xi)3,  {Xj)4,  {Xj)s,  (X!)s(  {x2)1(  (x2)2,  (x2)3,  (x2)4,  (x2)5,  (Xa ),-»>];  *) 
n 

eigenAve  =  f«  eigenewa  ; 

a=l 

Do[A[[i,  y]]  = eigenAve [ [i] ] ,  {i,  6}]; 

Do  [Aprime  [  [i,  y]  ]  =  edistnew2  [  [i]  ]  ,  {i,  6}  ] ; 

Do [distrainl [  [ i,  y]  ]  =  edistnewi  [ [ i]  ] ,  {i,  6}]; 

Do  [distrain2  [  [i,  y]  ]  =  cdistnew2  [  [ i]  ] ,  {i,  6}]; 

}.  (Y*  6}]; 

{*  Cnew=Cref  .W.  ( InvW+  (Sv-invW)  *W .A)  .  Inverse  [  InvW+Sv  .W.  A]  .InvW;  *) 

Cnew  —  Cref  •  W  •  (InvW  +  { Sy  —  InvW)  *  W  •  A  +  Dcracfc  •  W  •  Cn0crack  )  • 

Inverse  [  InvW  +  Sv  •  W  ♦  A  +  Dcrack  •  W  •  Cnocrack  ]  •  InvW; 

Aprimenew  =  Aprime ; 

Cref  =  (Cnew  +  Transpose  [ Cnew  ]  )  /2,* 

Cref  =  Cnew ; 

U  new  =  Cref  [  [4  ,  4]]; 

_  Cref  [[lx  2]] 

Vne”  "  2  (Cre£[[l,  2]]  +/Jne„)  ' 

p.  _  2  fJ-Tiew  ( 1  +  Vnew  )  „ 

Knew=  3(1-2  Vnew  )  # 

Clear [StrainDisturb] ; 

Sv  = MakeTwo [ PopEshSphere [ vnew ] ] ?  (*  Eshelby  tensor  of  outer  ellipsoid  *} 

Do [Sa  =  MakeTwo [PopEshSphere [vn] ]  ,  {a,  n-1}];  (*  Eshelby  tensor  of  micro-elements  *) 

{DcraCk2,  Ccrack}  =  PennyC  racks  tress  [Cnocrack  ,  vn/  jun,  fcracki ,  f  crack2 ,  fcrack3  ]  ; 

J,  {Iterations}]; 


{*  Do[ 

{Cfinai  =  (Cref  .W.  (InvW+  (Sv-lnvW)  . W. A)  +Acrack)  . Inverse  [InvW+Sv  .W.A+A4 ]  .InvW; 
Cref  =Cfinal  /  A^new  — Cref  [[4/4]]; 

vnew  =  2  (c^fYfi^nlLw)  ;  Sv  =MakeTwo  I  PopEshSphere  [  vnew  ]  ]  }  ,  { 1 0  }  ]  ;  * ) 
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c  final  -  Cref  f 

C final  =  (Cref  •  W  .  (  InvW  +  {Sv  “  IllVW)  •  W  .  A  +  DCrack  •  W  •  )  +  0  CcracJc  )  • 

Inverse  [  InvW  +  Sv  .  W  .  A  +  Dcracfc  .  W  .  Cnocrack  ]  •  InvW; 


MUdde  -  1  -  f i  ea  ..  _  2~ Sw  [  [ 4 >  4]]  ; 

Mref  L  1  J  J 


Mu^ds 


1  +  fl 


1 


f^ref  -^1 


1 _ 

-2  Sw  [{4,4)  ] 


Print [NumberForm [ {^1  ,  f  1  t  ^new  r  ^  ^ 

Print  [Mat rixForm[Cref  ]  ]  ; 

Print  [MatrixForm  [Cnocrack  ]  ]  ; 
comp  =  InvW  .  Inverse  [Cnocrack  ]  .  InvW; 
comp 2  =  InvW  .  Inverse  [cfinal  ]  .  InvW; 
Print  [MatrixForm  [Cfinal  ]  ]  ; 


MUdde  /  A^new  / 


MUdds  }  ]  ]  # 


Print  [  {c£lllal  [  [4,  4]  ] ,  ( 1  +  ferae  k3  ~  ^y-  j  j]; 

Print[{comp[  [3,  3]  ] /comp2  [  [3,  3]  ] ,  (l  +  fcrack3 - — i-tj  }]; 


Print [MatrixForm [comp2 ] ] ; 


Clear [OtherEigenAverage] ; 

Print  [MatrixForm  [  {cnocxack)  .  Inverse  [InvW  +  DcraCk  .  W  .  cnocraCk  ]  •  InvW]  ] ; 
Print  [MatrixForm  [  (Cnoctaok )  .  Inverse  [InvW  +  DoracIc2  •  W  .  Cnoorack  ]  .  InvW]  ] 
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Print  [MatrixForm  [  ( f x  Ci  +  £2  C2 )  *  Inverse  [InvW  +  Dcrack  ♦  W  »  Cn0crack  ]  •  InvW]  ] ; 
Print [MatrixForm [A]  ]  ; 

Print [MatrixForm [distrainl] ] ; 

Print [MatrixForm [distrain2] ] ; 

Clear[ "Global'*"] ;  (*  clears  all  variables  from  memory  *} 


Printed  by  Mathematica  for  Students 


